An itinerant half-metal spin-density-wave state on the hexagonal lattice 
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We consider electrons on a honeycomb or triangular lattice doped to the saddle point of the 
bandstructure. We assume system parameters are such that spin density wave (SDW) order emerges 
below a temperature Tjv and investigate the nature of the SDW phase. We argue that at T < Tn 
the system develops a uniaxial SDW phase whose ordering pattern breaks 0(3) x Z4 symmetry and 
corresponds to an eight site unit cell with non-uniform spin moments on different sites. This state 
is a half-metal - it preserves full original Fermi surface, but has gapless charged excitations in one 
spin branch only. It allows for electrical control of spin currents and is desirable for nano-science. 



Introduction: The electronic properties of single layer 
graphene have been the subject of considerable experi- 
mental and theoretical interest [H . Near half- filling, a de- 
scription in terms of non-interacting Dirac electrons cap- 
tures the essential physics, since interactions effects are 
suppressed by the low density of states (DOS). A sharply 
different behavior arises when graphene is strongly doped 
to 3/8 or 5/8 filling [2]. At this filling, a divergent den- 
sity of states and nested Fermi surface (FS) conspire to 
produce weak coupling instabilities to an extensive buf- 
fet of ordered states, including spin density waves (SDW) 
0) 01 , Pomeranchuk metals Q , and d wave supercon- 
ductors (SC) [i, i, i]. A similar situation arises on a 
triangular lattice at 3/4 filling [l^, U [ ■ 




It has recently been established using renormalization 
group (RG) methods @ that the two most relevant insta- 
bilities at weak coupling are towards SDW and a d-wave 
SC. The SDW vertex is the largest at intermediate RG 
scales, but superconducting vertex eventually overshoots 
it both at perfect nesting and away from perfect nesting, 
making d-wave superconductivity the leading weak cou- 
pling instability. The SC state has a d + id gap structure 
and breaks time- reversal symmetry [Hj . 

In this paper we assume that superconductivity is de- 
stroyed by an applied magnetic field, or alternatively that 
system parameters are such that corrections to the RG 
fiow become relevant before SC vertex overshoots SDW 
vertex. In both cases, the SDW instability becomes dom- 
inant, and an SDW order emerges. Previous work argued 
that the SDW state is non-coplanar and has a non-zero 
spin chirality 0, 0]. Such a state gaps out the entire 
Fermi surface (FS), i.e., is an insulator. The chiral SDW 
state has also been found in the strong coupling analysis 
for classical spins of fixed length [l^] • 

We argue that the situation is more complex than orig- 
inally thought, and the chiral SDW state is present only 
at the lowest temperatures. Over a wide intermediate 
range of temperatures, a different SDW state emerges 
in which SDW order develops simultaneously at three 
inequivalent wavevectors Qi, but the three vector order 
parameters are all aligned along the same axis. This state 
has a eight site unit cell with non-uniform spin moments 



FIG. 1: (Color online) a) The Fermi surface at the doping level 
of interest is a hexagon inscribed within a hexagonal Brillouin 
zone (BZ), for both honeycomb and triangular lattices. The 
FS has three inequivalent corners, which are saddle points 
of the dispersion, marked by a vanishing Fermi velocity and 
a divergent density of states. The three inequivalent saddle 
points Mi are connected by three inequivalent nesting vectors 
Qi, each of which is equal to half a reciprocal lattice vector, 
such that Qi = — Qi. b) Spin structure for the uniaxial SDW 
state. The SDW order quadruples the unit cell to a unit cell 
with eight sites (shaded). The enlarged unit cell has a large 
spin moment 3 A on two sites and a small spin moment —A 
on the other six. The total spin on each unit cell is zero. 



and zero net magnetization (Fig. [T|d). Such a state can- 
not be accessed starting from a spin Hamiltonian for local 
moments with a fixed length, and can only be accessed 
starting from a model of itinerant fermions. We show 
that in this state, unlike in any other known SDW state, 
the chemical potential shifts proportionally to the SDW 
order parameter preserving the original Fermi surface for 
one spin branch and gapping out the other spin branch. 
The uniaxial SDW state is therefore a 'half-metal' that 
allows for electrical control of spin currents. Such a state 
is highly desirable for nano-science applications. 

The model: For definiteness we focus on doped 
graphene at 3/8 filling. Our point of departure is the 
tight binding model [12], with the nearest-neighbor dis- 
persion 
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where fi — —ti at 3/8 filling. The FS then forms a per- 
fect hexagon inscribed within a hexagonal BZ (Fig. [1^). 
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The perfect nesting of the FS in doped graphene is quite 
robust - it is broken only by third and higher neigh- 
bor hoppings, which are generally quite small. The 
Fermi velocity vanishes near the hexagon corners Mi = 
(2^3,0), M2 = (7r/3,VV3), M3 = (tt/S, -tt/VS), 
which are saddle points of the dispersion: 

3ti 3ti 
^Mi+k ~ —(kySkl), eM2,3+k ~ — ^2fcj,(fcy=FV3fcj;), 

(2) 

where each time k denotes the deviation from a saddle 
point. Saddle points give rise to a logarithmic singularity 
in the DOS and control the SDW instability at weak cou- 
pling. There are three in-equivalent nesting vectors con- 
necting in-equivalent pairs of saddle points (see Fig.[T^): 

Qi = (0,27r/V3), Q2,3 = (±7r/3,-7r/\/3). (3) 

Each Qi is equivalent to — Qi modulo a reciprocal lattice 
vector. 

There are two electron-electron interactions that con- 
tribute to the SDW channel. One is a forward scattering 
interaction |k, k -I- Q^) — >• |k, k -I- Qi), while the other is 
an umklapp interaction, |k, k') — |k -I- Qi, k' -I- Qi). We 
label these interactions 172 and 173 respectively, for con- 
sistency with the notation introduced in Q. The parti- 
tion function for (72 — 93 model can then be written as 

Z = /i:'[V'^V]exp(-5[V'^V']), where S = J^^^ C(k,T) 
and 

a 

- .93V'l,a'0l,/3'0f.,/3'0f.,a " .92 V'a.a , (4) 

where the action is written in terms of electron operators, 
a, b are patch labels, and a and (3 are spin components. 
Each nesting vector has associated with it an SDW 



order parameter Aj 



The condition for the emergence of each is the same: 
((52 +ff3)Ai) log^ ti/TN = 0(1) [1, leaving a large num- 
ber of SDW states as potential candidates. We study 



the selection of the SDW order within Ginzburg-Landau 
theory and by comparing different SDW solutions in 
the mean-field approximation for Eq. (j4]) at arbitrary 
T<Tn. 

Ginzburg-Landau theory: To construct the Ginzburg- 
Landau theory, we decouple the quartic interaction terms 
by restricting the interaction to the spin channel and per- 
forming a Hubbard Stratonovich transformation to in- 
troduce the order parameters A^. We integrate out the 
fermions in the Matsubara frequency representation and 
obtain an action in terms of the order parameter fields 
Ai, which takes the form 
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For T fa Tn, we can expand ^ in small Ai/T/y. It is 
useful to define the expansion coefficients 
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where the integrands are expressed in terms of 
fcrmionic Green functions G — (iw„ — £k ^ Z^)^^? = 
(iw„-£k+Q. -^)"\ and G,,+j = (iw„ -£k+Q,+Q, -A^)"^ 
as 
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Ci = G^G 

= GG3G1G1+3, 



^2 — G^GsGi, 

C4 = G^GjGl (7) 



Diagrammatically, Zi -Z3 are given by 'square' diagrams 
with four ferniionic propagators and (Tap in the ver- 
tices, and Z4 is given by a 'hegagonal' diagram with six 
fermionic propagators, (see Fig. [2]) . The free energy eval- 
uated at T « Ttv can be expressed in terms of these 
coefficients as 



£ cx a(T - Tn) E + + ^2 + ^3)' + 2(^2 ~ Zi - Z^)[^\/^l + A^A^ + A^A^) 

i 

+ 4Z3((Ai • A2)2 -f (A2 • A3)2 + (A3 • Ai)2) - 4Z4(Ai • A2 x A3)' + • • • (8) 

I 



where a is an inessential positive constant. 

The quadratic term and the first quartic term in (jS]) 
set the overall magnitude of = X^i'^i; but do not 
differentiate between different SDW states. The second 
quartic term in ([5]) determines whether SDW order devel- 



ops only at one nesting vector, or at all three (depending 
on the sign of ^2 — .^1 — .Z3). Finally, the third quartic 
term and sixth order term control the relative orientation 
of the vector order parameters, if SDW order develops at 
multiple wavevectors. Close to the expansion to or- 
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ery site. Our result differs from earlier mean-field analy- 
sis which found non-coplanar insulating SDW order 



FIG. 2: (Color online) The terms quartic in A are pro- 
duced by processes represented diagrammatically by square 
diagrams. The diagrams for Z2 and Zz correspond to pat- 
terns A3,A3,Ai,Ai and A3,Ai,A3,Ai, respectively. The 
sixth order chirality sensitive term is produced by 'hexagonal 
diagrams.' Sample square and hexagonal diagrams are shown 
above. The integrals are dominated by momenta that bring 
all the fermion propagators to the vicinity of one of the saddle 
points of the dispersion. 



der is generally sufficient, but we include the sixth 
order term because .Z3 is suppressed by an extra factor 
of T/v/ii, which is exponentially small in the weak cou- 
pling limit. The relative smallness of arises because 
in the integrals for ^1,^2, and ^4, all fermions can be 
simultaneously brought to the saddle points, whereas in 
the integral for Z3, three fermions can be brought simul- 
taneously to saddle points, but the remaining fermion 
stays far away from not only the saddle point but also 
the FS. 

We evaluate the coefficients Z\ - Z4 to leading order 
in small T/v / 1\ and obtain [l5| 



Zx 



0.20 log ^ 
0.08 Tn 



0.58 
0.1 



(9) 



The positivity of Zi guarantees a second order phase 
transition, with the type of SDW order depending on the 
signs and relative magnitudes of Z2, Z3, and Z4. Since 
Z3 is smaller by T^/ti than Zi_2, and Z2 is smaller by 
log ^ than Zi, it follows that Z2 — Zi — Z^ < 0, so 
the system forms SDW order simultaneously at all three 
nesting vectors (the 3Q state). Meanwhile, the relative 
orientation of the three SDW order parameters is con- 
trolled by the sign of Z3 at the smallest A^, and by the 
sign of Zi at somewhat larger A^. Both Z3 and Z4 are 
negative and favor the non-chiral SDW order with the 
three A^ all aligned along the same axis. 

An order parameter of the form A(e**^^ '" + e^^^ '^ ± 
gjQ2 r^ leads to spin moments on the lattice of the form 
shown in Fig. [TJ A quarter of lattice sites have spin mo- 
ment 3A, the other three quarters have moment —A. 
Such an order cannot be obtained from any spin Hamil- 
tonian for local moments of constant magnitude on ev- 



at weak coupling. We note, however, the 3Q state that 
we found, with non-equal spin length on different sites, 
was not considered in that work and other earlier consid- 
erations of the type of SDW order. We found analogous 
results for fermions on a triangular lattice at Van Hove 
filling. This system is identical to graphene, except that 
the nesting is less robust and is spoiled already by second 
neighbor hopping. 

Properties of a uniaxial SDW: Is the uniaxial SDW 
state a metal or an insulator? To address this issue we 
need to compute the fermionic spectrum. Without loss 
of generality, we take the SDW to be uniaxial along the 
z axis, so that is a good quantum number, and spin- 
up and spin-down fermions decouple. Consider the state 
with Ai — A2 = A3 — AztTa. The up spins near the 
three Van Hove points are described by a simple 3x3 
Hamiltonian 
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e2,k - Sfj. A 
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where ei,e2,e3 are the dispersions near the Van Hove 
points, Eq. ([2]), and S^j. is the SDW-induced shift of the 
chemical potential. The 3x3 Hamiltonian describing the 
spin down branch is obtained by taking A — )• — A. At 
k = (i.e., at Van Hove points) the energies of spin- 
up excitations — Sfi are —A, —A, and 2A, and the 
energies of spin-down excitations are A, A, and — 2A. 
In conventional SDW states (e.g., SDW on a 2D square 
lattice) 5/x/A (X T^/Ef is negligibly small and can be 
safely neglected. We find, however, that in our case Sfi = 
—A, so that gapless excitations arise in the spin-down 
spectrum. 

To see the unexpected shift of the chemical potential, 
we diagonalize Eq. (jlOp and the corresponding equation 
for down spins and inspect six branches of excitations. 
We find that fixing (5/i — —A ensures that both in the 
paramagnetic and in the 3Q uniaxial SDW state there 
are four bands with £^k < and two bands with E^ > At 
for all momenta in the reduced BZ (see Fig. [3]). Since the 
chemical potential is fixed by the constraint that the total 
number of electrons (equal to the number of states below 
the chemical potential) must not change between A = 
and A ^ [l3|, it follows that we must set Sfj, = — A. For 
verification, we computed the thermodynamic potential 
fi(A,/z) from (O, numerically solved the simultaneous 
equations dfl/dA — and d^/dfi = ~N, and confirmed 
that Sfi — — A to a high accuracy. 

Having determined that Sfi — —A, we find from (jlO[) 
that gapless excitations emerge when £i,k£2,k£3,k = 0, 
which has solutions along three lines passing through 
each Van Hove point. Two of then coincide with the 
original FS, the third is directed towards the center of 
the BZ. The 3Q uniaxial SDW state is then obviously a 
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FIG. 3: (Color online) a) Excitation spectrum Sk ~ Ek — Sjj. of 
the 3Q uniaxial state. Negative k are along the FS, positive k 
are along the BZ boundary in the original BZ (along kx in the 
reduced zone). Placing the chemical potential at Sjj, = —A 
ensures that four bands lie below the chemical potential (hor- 
izontal dotted line) and two lie above for all k, irrespective of 
the value of A. Thus the choice n = —A conserves electron 
number. Excitations with spin projection opposite to A are 
in blue (solid), along A are in red (dashed) lines. Note that 
gapless excitations arise in the spin-down branch only, b) Free 
energy difference 5F = ^uniaxial — J'chirai between the 3Q uni- 
axial SDW state and the chiral state, evaluated in the mean 
field approximation for the honeycomb lattice Hubbard model 
with g2 = g3 = U = 1.7ti {Tn ~ 0.002ti). The 3Q uniaxial 
state has lower Free energy over a wide range of intermediate 
temperatures, but at the smallest T the non-coplanar, chiral 
state, studied in earlier works has lower Free energy. 



metal. We emphasize, however, that gapless states exist 
only for the electrons with spin projection opposite to 
A. The electrons with spin projection along A are fully 
gapped. Since a Fermi surface exists for one spin pro- 
jection only, we dub this state a 'half metal.' We found 
an analogous 'half-metal' spectrum for the 3Q uniaxial 
SDW phase on the triangular lattice. 

The half-metallic nature of the SDW should manifest 
itself in numerous experiments. For example, in tunnel- 
ing experiments conducted with electrons spin polarized 
along the z axis, a hard gap will be seen for down spins, 
but a Fermi surface will be seen for up spins. Further- 
more, since the low energy charged excitations involve 
up spins only, any charge currents will necessarily also 
be spin currents. Thus, the half metal state allows for 
electrical control of spin currents, which may be benefi- 
cial for nanoscience applications. 

Order parameter manifold: The uniaxial SDW order 
obviously breaks 0(3) spin-rotational symmetry. It also 
breaks discrete symmetry associated with either par- 
allel or antiparallel ordering of A^, i.e., in addition to the 
(A, A, A) state which we considered above, there are also 
(A, -A, -A), (-A,A, -A), and (-A,-A,A) states. 
These states have an identical structure of fermionic ex- 



citations, and correspond to the four in-equivalent ways 
to choose which two of the eight sites of the SDW unit 
cell carry large spins (see Fig. (Ub)). Equivalently, the 
three other states from the manifold are obtained from 
the (A, A, A) state by shifting the origin of coordinates 
to the center of one of three neighboring hexagons. An 
interesting possibility, which deserves further study, is 
that Zi symmetry can be broken before 0(3) symmetry, 
leading to a nematic-like state 

The phase diagram: Thus far we have constructed the 
Ginzburg-Landau expansion in small A/Tn- This expan- 
sion becomes less justified as we move towards zero tem- 
perature. To investigate the behavior at arbitrary T we 
calculate numerically the full Free energies of the various 
SDW states from ([5])- Upon doing this, we find that the 
3Q uniaxial state has the lowest Free energy over a wide 
range of intermediate temperatures, but undergoes a first 
order transition at a lower temperature to the insulating 
chiral SDW state discussed in earlier works iill^. We 
show the Free energy profile in Fig. [31d. We found this 
behavior both for graphene and for fermions on a trian- 
gular lattice. Intuitively, the chiral SDW state wins at 
the lowest T because it has spin-degenerate excitations 
and opens a full spectral gap, unlike the half-metal state. 

The Free energy profile in Fig. |31d is for weak/moderate 
coupling, when Tjv/ii ^ 1. At Tjv ~ ^i, the phase dia- 
gram is more complex and less universal. For complete- 
ness, we discuss the forms of Zi and the phase diagram 
at T/v ti in the supplementary material fl5|. 

Conclusion: We considered in this work the SDW in- 
stability on the honeycomb and triangular lattices, when 
doped to the saddle points of the dispersion. The SDW 
instability is subleading to a d— wave superconducting in- 
stability at weak coupling, but becomes the leading insta- 
bility if superconductivity is suppressed. We found that 
if the SDW ordering temperature is much smaller 
than the fermionic bandwidth, then a uniaxial SDW or- 
der develops simultaneously at three inequivalent nesting 
vectors. This has an order parameter manifold 0(3) x Z4 
and corresponds to the ordering pattern shown in Fig|T] 
Such a state can only be obtained from an model of itin- 
erant electrons with interactions, and not from a spin 
model of local moments. We found that such SDW state 
is a half-metal in which gapless excitations exist in one 
spin branch only. Such a state may be beneficial for 
nanoscience applications particularly because charge cur- 
rents will necessarily also be spin currents, which allows 
for electrical control of the latter. 

We thank L. Levitov for numerous discussions concern- 
ing the interplay between superconductivity and SDW 
order. We are also thankful to C. Batista, R. Fernan- 
des, I. Martin, and Fa Wang for useful conversations. 
G.W.C. is supported by ICAM and NSF-DMR-0844115, 
and A.V.C. is supported by NSF-DMR-0906953. 
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SUPPLEMENT 

In this supplement we present the calculations that were quoted in the main text. 

Calculation of Z\ 

We wish to evaluate 

^1 =^E/ ^G2(k,a;„)G2(k + Q3,c.„) (11) 

The integral over the Brillouin zone is dominated by those vahics of k where both Green functions correspond to 
states near a saddle point. Expanding the energy about the saddle points, we rewrite the integral as 

« TV / ^ ^ 2 (12) 

Where the integral is understood to have a UV cutoff for k of order 1. We now define a = ^/3tl/4{ky — V^kx) and 
b = ■^/Sti/'l{ky + V^kx), and rewrite the above integral as 

Z,=tV^/^^ ^ , (13) 

SVSh J-VtT (27r)2 (i^„ + ab)\ioJn - a{a + 6))' 

We now define x = ab and rewrite the integral as 

Z -- - ^ (14) 

We now assume Tjv <C t\ (which should certainly be the case for weak/moderate coupling). In this limit, we can 
perform the integral over x approximately, using the Cauchy integral formula, to get 

y _tV 2 /"^ rfg 1 2isigna;„ _ r^s;^ 4 rfg 1 isign^„(fl^ + 2ia;„)3 

The imaginary part of the above integral is odd in w and hence vanishes upon performing the Matsubara sum to leave 
an integral that is purely real 

_ sr^ 8|c.„| da 1 4cjl-3a^ 8|^„| da 1 Aujj 

' tr3V3tJ-V*r2^l«l(a^+«3 ~ ^3V3iii_^2^|a|(a4 + 4a;2)3 ^ ' 

with logarithmic accuracy. Performing the integral over a (again with logarithmic accuracy) gives 

El 1 ii I , U ^ 16.8 In Tp!- 
^^ln^ = ^ ' ^ 16.8 In + 10.5 (17) 

Where we take a;„ = 27r(n + l/2)Tjv, T = Tjv and perform the discrete sum on mathematica. 

Calculation of Z2 

This time we want to evaluate 

Z2=TY^J ^G\k,UJn)G(k + Q3,^n)G{k + Q^,Un) (18) 
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Again, wc anticipate this integral will be dominated by regions of the Brillouin zone where all three Green functions 
correspond to states near saddle points. Expanding the dispersion about the saddle points, we obtain 

z^^ry [ 4^ ^ i (19) 

^ J (27r)2 (i^^ _ 3|i (3^2 _ kl)y{iu:n - ^2ky{ky - V3k,)) (ioj^ - ^2ky{ky + V3k,)) 
Making the same coordinate substitutions as in the preceding section, we recast this as 

^ ^ (20) 

3^/3^l J-^ (27r)2 (^^^ + ab) (icjn - a{a + b)) (itOn - b{a + 6)) 
After scaling out Um, we can rewrite it as 

Z2=Ty — ^ / - — ^ =Ty ^ 2.9 (21) 

"^SVStilcJnl^ J -^{27t)^ (t + abY{i-a.{o. + b)){i-b{a + b)) ^ 12TT^V3ti\ujn\^ 

Where the rescaled integral is fully convergent, and can be done numerically on mathematica. The sum over Matsubara 
frequencies can also be done on mathematica, and yields the answer 

2.9x16.8 

Comparing with the previous expression for Zi, we see that Z2 ^ Zi x 2.9/ ln(T/ti). Thus, Z2 -C Zi provided 
the log is large. (If the log is not large then the evaluation of Zi with logarithmic accuracy does not suffice, and 
sub-logarithmic contributions to Zi must also be taken into account.) 

Calculation of Z3 

We want to evaluate 

Z3=T^J ^^G'(k,a;„)G(k + Q3,a;„)G(k + Qi,a;„)G(k + Qi + Q3,a;„) (23) 

This time it is not possible to place all the Green functions at the saddle points. In fact, we cannot even place all 
the Green functions at the Fermi surface - the best that can be done is to place three of the Green functions near a 
saddle point, but the fourth has to be off Fermi surface. Thus, we obtain. 



^s'^T^jj^-^— ITWT- M.^, /, T^TTTT- mTTTTT, : — T^TTrT- ^TT (24) 



1 

r)2 {icjn - ^(3A;2 - k^)) (iw„ - ^2ky{ky - T^k^)) {iun - ^2ky{ky + y^k^)){iLU - 2ti) 



Making the usual substitutions, and assuming ti Tjv, we obtain 

« V /*' (25) 

^ ^ l2s/ZT:HlLolJ-t^ {iuJn + ab){iuJn-a{a + b)){iuJn-b{a + b)) 

12VS7rHlLol J-t,/u,^ {i + ab) {i - a{a + b)) {i - b{a + b)) ^ ' 

The integral is convergent. As usual, the imaginary part is odd in to and vanishes and we care only about the real 
part. Performing the integral on mathematica and taking the real part, we obtain 

^ ^ V 6.5T;v ^ 6.5 

Which is parametrically smaller than Zi and Z2 by Tjv /ti . 
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FIG. 4; (Color online) The behavior of scaling functions fi{x) = fi (^77" j- ^) The scaling functions /i (black) and /2 (blue). 

b) The scaling function fs. c) The scaling function /s corresponding to the term Z2 — Zi — Z3 (solid line). Superimposed on 
this is a discrete plot of /i (points). Note that the scaling functions /a (a;), fiix) and fsix) change sign between small Tjv/ti 
(weak/moderate coupling) and Tjv ^ ti (strong coupling). 



Calculation of Z4 

We now calculate the coefficient of the sixth order chirality sensitive term in the free energy, u(Ai • (A2 x As))^. 
After some analysis of diagrams we find that, 

Zi^Tj2 /7^G2(k,^„)G2(k + Q3,L.„)G'(k + Qi,t^„) (28) 

We can now place all the Green functions on the near a saddle point. Making the usual substitutions, we obtain 

Z, = Y r^^'^" '^-^ (29) 

f;^ 12V37r2ti|w„|5 {i ^ ahf {i - a{a + h)f {i ~ b{a + hf) 

12V37r2ii|w,„|5 i_t,/T„ {i + ahf{i-a{a + h)f{i~b{a + hY) 

where in the second line we have assumed that the Matsubara sum is controlled by the first few Matsubara frequencies. 
The integral can be done numerically, and it is negative for small Tpf/ti. It then follows that at weak/moderate 
coupling Z4 < 0, so that the Free energy at sixth order also disfavors chirality. 



Scaling functions 

The above calculations were performed at weak/moderate coupling, assuming T^/ti <C 1. However, the inte- 
grals can be evaluated at arbitrary T^/ti. To this end, it is useful to define the scaling functions /{(T^/ti) — 
Zi{TN /ti) /Zi{0). These scaling functions are evaluated numerically and shown in FigHl 



PHASE DIAGRAM AT STRONG COUPLING 



At strong coupling, when (72,3 > ti and Tn ~ ti, our analysis based on Ginzburg-Landau expansion is less accurate 
because fermions can no longer be approximated as free particles (the self-energy corrections to fermionic lines and 
vertex corrections to square and hexagonal diagrams are generally of order one). Nevertheless, if we apply our analysis 
to Tff ~ ii, we find that Z3, Z4 and Zi — Z\ — Z3 all change signs at some Tm/H (see Fig. 

The first sign change occurs in the sixth order chirality sensitive term Z4, which becomes positive for /t > 0.1. 
When is positive, the chiral SDW state Q is energetically favored, provided we are sufficiently far below for 
the sixth order term to dominate over the quartic term Z^. Thus, at large Tjv/ii, the uniaxial SDW phase has a much 
narrower region of stability, and the transition into the chiral SDW phase happens quite close to T/v. 

The term Z\ — Z2 — Zy, is next to change sign, becoming negative for T^/T > 0.35. Once this term becomes 
negative, the system prefers instead a IQ collinear state, of the form discussed in [ll|, wherein SDW order develops 
only at a single nesting vector. The subsequent sign change of Z^ at T^q jt k, 0.55 has no physical consequences. 
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The IQ collinear SDW state that forms at T/v/t > 0.35 is a (full) metal because the entire FS is not gapped out. 
The competition between a metallic collinear state and non-coplanar insulating state has been detected numerically in 
the mean- field analysis at strong coupling , and our results for the strong coupling case are in line with this earlier 
study. Our strong coupling results are also consistent with the studies that found a non-coplanar, chiral SDW order 
in the models of spins of the same fixed length at every lattice site ^ [lO||. However, in the weak/moderate coupling 
limit, our results indicate that the preferred state is a uniaxial SQ state of a sort not considered before, which can 
only be realized in a model starting from itinerant fermions. 

We also analyzed the evolution of Zi with T/v/ii for fermions on a triangular lattice. We found similar trends, e.g., 
sign change of Z4. However, for a triangular lattice, the first sign change (in Z4) occurs at a much larger Tpf/ti ^ 0.5, 
when the itinerant approach is very questionable. 
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